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We consider physical properties of a superconductor with a recently proposed type of odd- 
frequency pairing that exhibits diamagnetic Meissner response (“odd-dia state”). Such a state was 
suggested in order to address stability issues arising in an odd-frequency superconducting state with 
paramagnetic Meissner response (“odd-para state”). Assuming the existence of an odd-dia state 
(due to a proper retarded interaction), we study its coexistence with an odd-para state. The latter 
is known to be generated as an induced superconducting component in, e.g., singlet superconduc- 
tor/ferromagnet proximity structures or triplet superconductor/normal metal systems. Calculating 
the superfluid density of the mixed odd-para/odd-dia state and the Josephson current between 
the odd-para and odd-dia states, we find that the expressions for the currents in both cases have 
non-vanishing imaginary contributions and are therefore unphysical. We show that a realization of 
the odd-dia state implies the absence of a Hamiltonian description of the system, and suggest that 
there exists no physically realizable perturbation that could give rise to the spontaneous symmetry 
breaking necessary for an actual realization of the odd-dia superconducting state. 
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I. INTRODUCTION 

Pairing of fermions in superconductors and superfluids 
can be described by the so-called anomalous Green func¬ 
tion F(l;2) (here 1 and 2 denote sets of coordinates for 
the two fermions). The Pauli principle requires antisym¬ 
metry under permutation of two fermions in a Gooper 
pair, 1) = —F(l; 2), leading to the standard classifi¬ 
cation [l| of superconducting phases: even parity under 
exchange of the spatial coordinates of F (s-wave, d-wave, 
etc.) must be accompanied by odd parity under exchange 
of the spin coordinates (singlet). Similarly, odd parity 
in the spatial coordinates (p-wave, /-wave, etc.) must 
be accompanied by even parity in the spin coordinates 
(triplet). This classification scheme assumes that F stays 
unchanged under permutation of two imaginary-time co¬ 
ordinates (we assume for definiteness the Matsubara rep¬ 
resentation) . 

In 1974, Berezinskii pointed out Q that this classifi¬ 
cation, implicitly assuming even dependence of F on the 
imaginary time difference t = ti — T2, can be doubled if 
one allows for an odd symmetry with respect to time (or 
to the Matsubara frequency ui in the Fourier represen¬ 
tation). This opens up a possibility of, e.g., an s-wave 
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triplet ii or a p-wave singlet Q phase without vio¬ 
lating the Pauli principle or the Fermi-Dirac statistics. 
These exotic phases are called odd-frequency (odd-w) 
states. Although there has been much theoretical work 
on odd-frequency pairing states (as the principal state, 
i.e., the state corresponding to a normal-state instab ility 
in the pairing interaction channel) in bulk systems |5I-I2(tI| , 
experimental evidence is still lacking. 

However, it turns out that the s-wave triplet state can 
be realized as an induced (as opposed to spontaneously) 
symmetry-broken state, leading to an induced supercon¬ 
ducting component in a rather conventional system con¬ 
sisting of an s-wave singlet superconductor (S) and a 
ferromagnet (F) [l^. Ferromagnetism here breaks the 
symmetry in the spin space and opens up the possibility 
of generatin g triple t superconducting correlations from 
singlet ones [271 - 14^ . Alternatively, if the triplet corre¬ 
lations are already present in the system due to, e.g., a 
principal even-w p-wave triplet state (TS), then an odd-w 
s-wave triplet state can be generated due to breaking the 
isotropy in real space This can be achieved, 

e.g., by contacting the superconductor with a diffusive 
normal metal (DN) |^. Odd-frequency supercon¬ 
ductivity in such systems is an emergent inhomogeneous 
phase appearing in the vicinity of an interface and pene- 
tratin g in to the adjacent materials due to the proximity 
effect 0,113,111. 


Induced odd-frequency correlations have also been dis¬ 
cussed in connection with surface Andreev bound states 
in unconventional superconductors (lol - IH^ . vortex sys¬ 
tems [SSI - ISTII. Mai orana fermions in topolo gica l super¬ 
conductors [SSHg^, and multi-band systems 651. 

An unexpected feature of the induced s-wave triplet 
state in diffusive proximity structures is that the resulting 
local density of state has a zero-energy peak [1^ |H-lz3 ■ 
For ballistic structures, where both even-w (e.g., p-wave) 
and odd-w (e.g., s-wave) triplet states are generated at 
the superconductor-ferromagnet interface, a low-energy 
Andreev bound state band around zero energy was found 
[iSlnl. Furthermore, the sign of the local current re¬ 
sponse to an external vector potential turns out to be 
unconventional (in other words, the local superfluid den¬ 
sity is formally negative, ns < 0) This would 

lead to an instability of such a principal superconduct¬ 
ing state in a bulk homogeneous system [34lfll. [illTii. 
On the other hand, the induced odd-w state is stabi¬ 
lized in inhomogeneous proximity systems due to pres¬ 
ence of a principal conventional superconductor (note at 
the same time that under special conditions instability in 
thin-film proximity structures is also possible [Z^). Since 
the conventional {ns > 0) superconducting response cor¬ 
responds to the diamagnetic Meissner effect, we refer to 
the unconventional {ns < 0) state as a (locally) paramag¬ 
netic one (odd-w-para state) SHI. At the same time, 
in order to avoid confusion, we note that if one studies 
the Meissner effect in a proximity structure where the 
surface region has a prevailing odd-w-para component 
, then this surface region would show an oscillat¬ 


ing Meissner effect [tI, [tJ (since the sign of ns enters 
the screening length under the square root) which finally 
turns to a diamagnetic one in the bulk of the conventional 
superconductor. 

Although surprising at first sight, the odd-w-para state 
is predicted on the basis of microscopic models and is un¬ 
doubtedly realized as an induced superconducting state 
in real proximity structures [23 - 1^ . There are numer¬ 
ous experimental results consistent with manifestations 
of this state . An interesting feature of the odd-w- 

para state in proximity structures is that it can be spa¬ 
tially separated from other superconducting components 
arising due to the proximity effect p3 - l43| . For example, 
in S/F junctions the so-called long-range proximity effect 
arises [0|, meaning that not too close to the interface 
the only surviving superconducting component in the F 
layer has odd-w-para properties. So, this superconduct¬ 
ing state is indeed realized (as the only superconducting 
component) in certain spatial regions. 

Recently, a further extension of the classification of su¬ 
perconducting phases was proposed in Refs. [13, [^. The 
proposed state is also odd-w, but at the same time it 
possesses a conventional diamagnetic Meissner response 
(odd-w-dia, ns > 0), and thus is supposed to solve the 
problem of instability for a bulk material. If so, this 
state could be realized as a homogeneous principal su¬ 
perconducting phase. Being more conventional at first 
glance, this proposal relies on an essential assumption 
that a proper retarded interaction leading to this state 
can exist in a real material. 

Motivated by this proposal, we study physical proper¬ 
ties of the odd-w-dia state (assuming it is actually real¬ 
ized). It turns out that while the new state seems con¬ 
sistent by itself, coexistence of the two different odd-w 
states (para and dia) leads to unphysical consequences. 
Below, we describe the inconsistencies arising from the 
assumption of such coexistence and propose our view on 
a possible resolution of the contradiction. 

The paper is organized as follows. In Sec.|Bl we review 
the path-integral approach of Refs. [ 13 , in order to es¬ 
tablish the notations and underline the features that are 
essential for further discussion. In Sec. IBI Al we consider 
possibility of coexistence of the odd-w-para and odd-w- 
dia states in the same spatial region. In Sec. IBIBl we 
consider the possibility of their coexistence in the Joseph- 
son junction. Our results are then discussed in Sec. IIVI 
and summarized in Sec. El 


II. PATH-INTEGRAL FORMULATION 

In this section, we develop the path-integral formalism 
for describing superconductivity in both homogeneous 
and inhomogeneous states with general interactions. We 
then proceed to discuss a specific type of interactions, 
following the path-integral approach of Refs. [93, for 
describing the homogeneous principal odd-w-dia state. 
Reference [93| treated spinless fermions, while the gen- 
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eral case of arbitrary pairing for fermions with spin was 
considered in the Appendix of Ref. . We briefly out¬ 
line the main points of the derivation, following Ref. [9^ , 
with intention to set up the general framework necessary 
for our later discussion of coexistence between the odd- 
w-dia and odd-cu-para states. We pay special attention 
to differences between description of the odd-w-dia and 
-para states. 

The reader less interested in the technical details of 
a rigorous treatment can skip the following subsection 
which offers a general treatment, and proceed directly to 
the next subsection that presents a simplified treatment 
valid for a more restricted set of interactions. 

We adopt the following definition for the complex con¬ 
jugation of Grassmann variables: '0** = '*/'> ( 0 i 02 )* = 
02 01 - The temperature is denoted as T, and we put 
h = /cb = 1 throughout the paper. 

A. General treatment 

1. Interaction 

We consider a general interaction 

V0/3;75(riri,r2T2;r3r3,r4T4) (1) 

as matrix R(12;34) with compound index 12 denoting 
the set of variables riTia, V 2 T 2 I 3 and compound index 34 
denoting the set of variables rsrsy, v^taS. Here, a, /3, 7 ,5 
are spin indices, while the imaginary times (we use the 
Matsubara representation) are denoted ri,...,r 4 , and 
the spatial coordinates are ri,...,r 4 . We use a brief 
notation for integration 

/ =//, etc. ( 2 ) 

Ji Jo J ^ J12 J1J2 

We will use the singular value decomposition for the in¬ 
teraction matrix defined as above. Accordingly, the inter¬ 
action matrix can be expanded into the singular vectors 

y(12; 34) = ^ C^(12) Vx 7^34)0 (3) 

where A labels a complete orthonormal set of left singular 
vectors C''( 12 ), 

[ C^(12)*C^'(12)=(5av, (4a) 

J12 

^C^(12)C\34)* = di3-524, (4b) 

as well as correspondingly for the right singular vectors 
77 ^( 34 ). We use a brief delta-functional notation, e.g., 
03 = <5(ri - r3)(5(ri - T3)Saj. 

The singular values V\ are real and can be assumed all 
non-positive (we slightly modify the convention for the 
usual definition in terms of non-negative singular values). 


We will also use the pseudo-inverse of this matrix, defined 
as 

R-P(12;34) = J 7^(12) C^(34)*, (5) 

where is obtained from V\ by replacing all non-zero 
singular values by their inverses, leaving all zero singular 
values untouched. 

Below, for brevity we employ the matrix product no¬ 
tation 

[AR](12;34)= [ A(12; 56)5(56; 34), (6a) 

J56 

[aA](34)= / a(12)A(12;34), (6b) 

J12 

[Aa](12)= / A(12;34)a(34), (6c) 

J34 

for dealing with four-point functions (A, B) and two- 
point functions (a). 

Then Pi (12; 34) = [R1 /-p]( 12; 34) and P2(12;34) = 
[1/-PR](12; 34) turn out to be Hermitian projec¬ 

tor matrices (not necessarily identity matrices), and 
[V-PVV-P] = v-p as well as [VV-pV] = R. 

2. Partition function and action 

The partition function Z is written in the path-integral 
formulation with the help of Grassmann fields 0(1) and 
001 ): 

Z = J P0*P0 exp(-5o - 5i„t), (7) 

5o = ^001)(9.,+O0(l), (8) 

Sint = l [ p012)R(12; 34)034), (9) 

^ dl234 

where ^ = —d^^l2m — /x is the kinetic energy counted 
from the chemical potential, and notations 

p(12)=0(l)0(2), p012)= 002)001) (10) 

are introduced for brevity for the pair-density field and 
its complex conjugate. 

The So action describes free particles, while 5int de¬ 
scribes spin-dependent interaction. Requiring Z = 

Z* and changing the integration variables (Grassmann 
fields) into r]{r,T,a) = 0(r, l/T — r, a), ?7*(r,r,a) = 
0*(r, 1/T — r, a), we obtain the condition 

R(12;34) = 0034;I2), (11) 

with 1 = (ri, IjT — ri,a), etc. This is equivalent to 

Vx f 0(12)*0'(I2) = 14' / C^12)*r?^'(l2). (12) 

J12 J12 
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At the same time, V fulfills the relations 


4- Saddle-point solution 


F(12;34) = y(21;43) = -^(21; 34) = -^(12; 43), 

(13) 

conditions that follow directly from exchanging integra¬ 
tion variables in iSint ■ These conditions are equivalent to 
the relations 

C^(12) =-C^(21), 77^(12) = V(21), A^Ao, (14) 

where Aq is the set of all A for which Va = 0. With the 
definitions 

px= [ r7^(12)V(12), p+= [ p*(12)C^12), (15) 

Jl2 Jl2 

we obtain finally 

= l f ptVx PA. (16) 

3. Mean-field approximation 


Now, we focus on finding the saddle-point solution. 
Instead of integrating over A and A, we take a trial path 
with respect to which we will vary the action. Then, we 
can define the mean-field free-energy functional: 

J'MpiA, A] = —TIuZmf 

= -Tln J D^p*D^j;exp{-SMFbP*,tp,^,^]) , (23) 

where 5 mf = Nq -f 5aux + Sa is the mean-field action 
in which A and A are replaced by the trial path. To 
find the saddle point, we should minimize J^mf with re¬ 
spect to the trial path. As a result, we obtain the self- 
consistency equations for the mean-field pair potential 
Amf (also called the gap function), 

Amf(12) = [PiA]( 12) = - / F(12;34)F(34), (24a) 

J 34 

Amf(34) = [AP2](34) = - / F+(21)y(12;34), (24b) 

J12 


Performing a Hubbard-Stratonovich transformation 
amounts to multiplying the partition function by a finite 
constant 


DAD A exp 


/ A(12)F"P(12;34)A(34) L (17) 

U234 J 


where we have introduced the anomalous averages 
(anomalous Green functions), 

F(12) = (p(12))j^p, F+(21) = (p*(12))^p, (25) 

with the help of the definition of the mean-field averaging: 


and subsequently shifting the integration fields as A >->• 
A -I- \yp] and A H> A -|- [p*V\, using [(A -|- p*V)V~P{A -\- 
Vp)] = [AtA-PA] -h [p*PiA] + [AP 2 P] + [p*Vp]. As a 
result, the four-fermionic term 5int is decoupled: 

exp(-5i„t) J DAD A exp(-5aux - Sa), (18) 


^ / g7 /.* iJ7 /;(...) e -^ MP [ V ^->, A , A ] 

The notations F and F'^ for the two anomalous averages 
are standard in the theory of superconductivity, and the 
-b superscript should not be confused with the Hermitian 
conjugation f. From the relations 


where 

5aux = -i^Jp*(12)[PiA](12) -b [AP2](12 )p(12)), 
Sa = -1 [ [AP2](12)TA-^’(12;34)[PiA](34). (19) 

Defining the fields 

Aa = / C^(12)*A(12), Ax= [ A(12)77^(12) (20) 
J12 J12 

for all A ^ Aq, and setting Aa = A a = 0 for all A G Aq, 
we obtain 

5aux = -^ / (p^Aa + Aapa), (21) 

SA = -ll^(AxVr^x)- (22) 

We require Aa = A^ for ensuring convergence in the 
Hubbard-Stratonovich transformation. 


i^(12) = -J^(21), i^+(12) = -i^+(21), (27) 

following directly from the properties of Grassmann vari¬ 
ables, and from Eqs. m, one obtains the corresponding 
relations 

Amf(12) = —Amf(21), Amf(12) = —Amf(21). (28) 

Finally, with the definitions 

Fa = / 77^(12)*F(12), F+ = [ F+i21)C\l2) (29) 
J12 J12 

(note the difference in the definitions for F\ and Aa , and 
for F^ and Aa, correspondingly), we obtain 

Aa = -VxFx, Al = -VxF+. (30) 

From these relations it follows that, for all components 
Fa, F^ for which A ^ Aq, the symmetry relation F^ = 
F^ holds (remember that Vx is real). For all Fa, F^ 
with A G Aq nothing follows from these equations (as 
then 14 = Aa = A^ = 0). 
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B. Simplified treatment 


spin indices (a, /3, 7 , d) is implied, and notations 


1. Interaction 


In order to make our consideration more transparent 
and facilitate comparison with Refs, we now as¬ 

sume that the interaction is homogeneous with respect 
to spatial coordinates and time. For that we slightly 
change notations, indicating the spin indices explicitly, 
while the spatial coordinates and imaginary times are 
gathered in 4-vectors denoted as 1 = Xi = (ri,Ti) and 
2 = X 2 = (i’ 2 ,t- 2 )- The general interaction of Sec. Ill Al 
[Eq. ([T])] now takes the form —2)(5(1 —3)<5(2 —4), 

where we denote i5(l — 3) = i5(ri — r 3 )( 5 (ri — T 3 ), etc. 

At the same time, we keep general form with respect 
to interactions in both the singlet and triplet channels 
(so that later we can consider specific cases on the basis 
of general equations). So, we assume interaction of the 
form 


Va,0,^sil-2) = Vsil-2) 




^y«(i-2) 

3 = 1 


{iaja2)ctp 


(31) 


The singlet interaction 14(1 — 2) and the three compo¬ 
nents of the triplet interaction {1 — 2 ) are assumed to 
be either negative definite or zero [Q^. For this case, the 
matrix inverse of the interaction matrix with two com¬ 
pound spin indices is defined as 


[V-\l-2)]^p,^s 


Vs{l-2)-^ 


{i(J2)ap{i(T2)*^S 

2 


-f ^E/^41-2)-1-—^ ^ (32) 

i=i 

with V{1 — 2)~^ = l/R(l — 2) for nonzero components 
and zero otherwise. 


2. Partition function and action 

The partition function Z is written in the path-integral 
formulation with the help of Grassmann fields ipai^) 

C(l): 

Z = J Sint), (33) 

5o = ^C(l)(9..+OV'c.(l), (34) 

5int = ” 2)p:^(1, 2)p^s{l, 2), (35) 

where ^ = —dn^j2m — p, is the kinetic energy counted 
from the chemical potential, summation over repeated 


P„;3(1,2)=V'„(1)V^;3(2), p:^(1,2)=V'^(2)V^:(1) (36) 

are introduced for brevity for the pair-density field and 
its complex conjugate. The brief notation for integration 
is 



The So action describes free particles, while 5int de¬ 
scribes spin-dependent interaction. Requiring Z = 
Z* and changing the integration variables (Grassmann 
fields) into r]a{x) = tpai—x), rjaix) = '0*(—x), we obtain 
the condition 

K./3;^5(1-2) = R4„^(2-1). (38) 

At the same time, V fulfills the relation 

Ra/3;75(l-2)=R^„;5^(2-l), (39) 

a condition that follows directly from exchanging inte¬ 
gration variables in 5int. These two symmetries lead to 
the requirements that the singlet and triplet interactions, 
14(1 — 2) and V'/'^^(l — 2), are real and even: 

K(l-2) = 14(2-1), r/^')(1-2) = r/^'^(2-1). (40) 

3. Mean-field approximation 

Performing a Hubbard-Stratonovich transformation 
amounts to multiplying the partition function by a finite 
constant 

X A7(l,2)A,j(l,2)|, (41) 

and subsequently shifting the integration field as 
Aa/3(1, 2) H> Aa/3(1, 2) -|- Vap-,-ys{l — 2)p^s{l, 2). 

As a result, the four-fermionic term Sint is decoupled, 

exp(-5int) J DA*DA exp(-5aux - 5 a), (42) 

with 

5aux = / [^«/3(l, 2)p„^(l, 2 )-|-Aq,^(1, 2)pq;)3(1, 2)] , 

J 12 

(43) 

5a = [R-'(l - 2)]^^^^, A;^(1,2)A^,(1,2). 

(44) 

Restricting the Hubbard-Stratonovich transformation to 
negative definite 14(1 — 2) means that the quadratic form 
in the exponent of Eq. m is negative definite, ensuring 
that the integration in Eq. m is convergent. 
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J^. Saddle-point solution 


Now, we focus on finding the saddle-point solution. 
Instead of integrating over A* and A, we take a trial path 
with respect to which we will vary the action. Then, we 
can define the mean-field free-energy functional: 


-?'mf[A*, A] = —TlnZMF 

= -rin J Dip*DipeTq>{-SMF[ip*,ip,^*,^]), (45) 

where iSmf = iSq -f 5aux + Sa is the mean-field action 
in which A* (1,2) and A(1,2) are replaced by the trial 
path. Further below we will treat the homogeneous case, 
for which the saddle-point solution does not depend on 
the center-of-mass coordinate (cci -I- X2)f‘2] for this case 
one chooses the trial path A*(a:), A(x) depending only 
on the relative coordinate x = Xi — X 2 - 

To find the saddle point, we should minimize Juvif with 
respect to the trial path. As a result, we obtain the self- 
consistency equations for the mean-field pair potential A 
(also called the gap function): 


A,^(l, 2) = - 2)F^5(1, 2), (46a) 

A:^(1,2) = - l)F+(2,l), (46b) 

where we have introduced the anomalous averages 
(anomalous Green functions), 

F^p{l,2) = (47a) 

F+^(1,2) = (C(1)V’^(2))mf> (47b) 

with the help of the definition of the mean-field averaging: 

The notations F and F+ for the two anomalous averages 
are standard in the theory of superconductivity, and the 
-|- superscript should not be confused with the Hermitian 
conjugation f. Grassmann variables ensure fermionic an¬ 
tisymmetry of the anomalous averages: 

F„^(1,2) = -F;3„(2,1), F+^(l,2) = -F+„(2,l). (49) 

Together with the symmetry of the interaction, Eq. (15^ , 
this property translates into fermionic antisymmetry of 
the pair potential, 

A„;3(1,2) = -A;3c(2,1). (50) 


The M matrix is written explicitly in the particle-hole 
space, and each its element is a matrix in the spin space. 

In addition to the Gor’kov (anomalous) Green func¬ 
tions (|T7|). we also introduce the standard Green func¬ 
tions: 


G„;5(l,2)=-{V^„(l)V^^(2))^p, (54) 

2) = -{rM)M^))^F = -Gpo.i2, 1). (55) 

The quadratic form of the action (15^ implies that the 
Green functions are expressed in terms of the M matrix 


Ma/3(1,2) 


(-Ga.pil.2) F„;3(1,2) \ ' 

1,F+^(1,2) -G'^^{h2)) 


(56) 


C. Symmetry relations 

In this subsection we discuss a relation between the 
two types of the anomalous averages, Eqs. (07]), for sim¬ 
plicity focusing mainly on the homogeneous case, when 
the Green functions depend on the difference of the co¬ 
ordinates, (I — 2). In the homogeneous case, Fourier- 
transformed functions A(l,2) [see Appendix lA] take the 
form A{k, k') = A{k)6 {[k — fc']/27r). 

Similarly to the procedure discussed in Ref. [93| (gen¬ 
eralizing it to the case of fermions with spin), we can 
obtain a relation between F and F+ directly from def¬ 
initions (I25|) and (El), applying complex conjugation 
to one of the anomalous averages. This procedure is 
nontrivial because in the path-integral formulation with 
averaging defined according to Eqs. dMl) and dm) , we 
have 5 mf due to the dr term in the action 

5mf, and then (A)]y[p is not necessarily equal to (A*)]y[p 
in the general case. In order to relate F* to F'^, 
we have to take into account fermionic antiperiodicity, 
ipaiT + /3) = —ipaiT), change r —r, and define new 

variables of the path integration rj depending on the sym¬ 
metry of A: for the case of A(t) = A(—r), we define 
r]air,T) = V^a(r,-r), ry*(r, r) = while for the 

case of A(t) = — A(—r), we define ??a(r, r) = 
r?*(r,T) = —z')/;* (r, —t) [q^- The relation then depends 
on the symmetry of A(t) or A(w). We find 

F+^(k,a;) = SAi";jk,-cu), (57) 


Introducing for notational simplicity 

A+^(l,2) = A;„(2,I), (51) 

we can write the fermionic part of the action in the form 


So 


Mc^{l,2) 


i^^(C(l),V’a(l))M„/3(l,2) 

^ (52) 

fSil-2)S^p{dr+0 A„^(1,2) \ 

A+^(1,2) S{l-2)S^p{dr-0)' 

(53) 


where sa = ±1 for the even-/odd-a; dependence of A. 
The type of Meissner response is eventually determined 
by the relative sign between F^p{k) and Fp^{k), so it 
depends on both sa and the frequency symmetry of the 
anomalous averages. 

In the model of Refs. [ 13 , [S^, it was implicitly assumed 
that the frequency symmetry of A directly determines 
(coincides with) the frequency symmetry of F and F^. 
In this case, Ftj/ 3 (k, —w) = SAA'a,g(k, w), and then Eq. 
(1571) immediately yields F^(fc) = Fp^{k) (both in the 
even-w and odd-w cases). This relation corresponds to 
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the conventional diamagnetic Meissner response, there¬ 
fore the authors of Refs, concluded that the odd- 

w superconducting state can be realized as a bulk state 
(if a proper interaction leading to odd-w pair potential 
really exists). 

On the other hand, relation (EZD is more general and 
remains valid, for example, if the Zeeman term (exchange 
energy) is added to the action determined by Eqs. (l34l) 
and (j35p . This term breaks the symmetry in the spin 
space and leads to appearance of superconducting com¬ 
ponents (anomalous averages) with the symmetry differ¬ 
ing from that of the interaction V and the pair potential 
A. The simplest example is the “superconducting ferro- 
magnet” considered in Appendix |B] Since the (instanta- 
neous) interaction in this model is of conventional BCS 
type [99| , the pair potential is defined at coinciding imag¬ 
inary times so that instead of A(r) we only have nonzero 
value A(0). This situation is only compatible with the 
symmetry class having even-r dependence, for which the 
case of sa = 1 is realized. However, the anomalous av¬ 
erages are mixtures of even-w singlet and odd-w triplet 
components [see Eqs. (lB4l)-(lB6l)]. So, while A is even-w, 
F and F~^ do not have definite frequenc y sy mmetry. 

Similar sit uati on is realized in S/F and TS/DN 

113, Il0l| proximity structures. Generally, this 
type of odd-w states does not correspond to the symme¬ 
try of the underlying effective electron-electron attrac¬ 
tion, and is induced due to symmetry breaking either in 
spin or coordinate space. At the same time, A belongs 
to the even-w class so that relation (EH with Sa = 1 is 
realized. This relation is valid for all components of the 
anomalous averages. For odd-w components, it yields 
~ which corresponds to the paramag¬ 

netic Meissner response. 

Note that relation (l57ll was obtained directly from the 
definitions of the anomalous averages. One may suggest 
that an alternative way to obtain a relation between F 
and E+ could rely on the self-consistency equations (|46|) . 
In order to avoid confusion, we note that in general the 
pair potential can be directly calculated from F and the 
interaction V, however, for the calculation of F in terms 
of A in general a set of differential equations must be 
solved (the Gor’kov equations), which depends on addi¬ 
tional interactions and, importantly, on boundary con¬ 
ditions. If one decomposes the functions F and F~^ into 
symmetry components, F = F+F, where F has the same 
symmetry as A, while F has differing symmetries [see, 
e.g., Eq. (Ull)], then A^/s = -Va/3;js(l - 2)Fjs(l,2) and 
Ra/ 3 ; 7 i 5 (l — 2)E^5(1, 2) = 0 hold. Then, only for the com¬ 
ponent F can one obtain a relation F^^(l, 2) = Fp^{2, 1), 
finally leading to E+^(fc) = F^^{k). This relation ap¬ 
plies, e.g., to conventional bulk superconducting states, 
or the odd-w-dia state of Refs. [93, However, it 

does not apply to the induced odd-w-para components of 
Refs. [ 23 ,0, [ 33 ,1131 , since about the relation between F 
and F^ nothing can be inferred from the self-consistency 
equations. 

We call the superconducting correlations with the same 


symmetry as the pair potential principal components, 
while other components — induced (e.g., by additional 
interactions, interfaces, other inhomogeneities, or exter¬ 
nal fields). In other words, the principal components of 
F and F~^ correspond to the symmetry components of A, 
while induced components of F and F'^ do not contribute 
to the right-hand sides of Eqs. (|151) due to the structure 
of V. The odd-w-dia state discussed in Refs. [13,1111 was 
implicitly assumed to be the principal component. 


D. Homogeneous case 

In the homogeneous case, the self-consistency equation 
(I46a|) takes the form 

A„/3(fc) = - Jidk')Vai3-^s{k - k')Fjs{k'). (58) 

Relation m turns to (IA91) and is simplified as 

A+^ik) = A}^{k). (59) 

In the general case, the pair potential can be decom¬ 
posed into the singlet component do(fc) and the triplet 
component d(fc) as 

A„/3(fc) = do (fc) ( 10 - 2 ) 0/3 + d(A:)(fcro-2)o/3, (60) 

where ct/ with i = 1,2,3 denotes the Pauli matrices in 
the spin space (below we will also use o-q to denote the 
unity matrix). Fermionic antisymmetry of A [Eq. (1501) ] 
implies 


do{k) = do{-k) = SAdoi-K^), (61) 

d(/c) = —d(—/c) = —SAd(—k, w). (62) 


When discussing principal odd-w states with sa = — 1 
(recall that superconducting components of A are prin¬ 
cipal by definition), we assume that it can indeed be re¬ 
alized due to some proper interaction. The possibility to 
find such an interaction is not at all clear to date. 

From Eq. (1551) we then find the Green functions [^ : 

Gap{k) = —G'p^{—k) = Go{k)Sa/3 + G(/c)(Tq,^, (63) 

Fapik) = F^^{k)* = Fo(fc)(fCT2)a/3 -f F(fc)(icro-2)a/3, 

(64) 


where the scalar and vector components are given by 


p \ (*w -I- ^k)[w^ + + F>oik)] 

[oj^ + El{k)][uj^ + Elik)]' 


G{k) = 
Fo{k) = 


F{k) = 


_(icu -I- Ck)D(fc)_ 

[uj^ + Elik)] [uj^ + El{k)] ’ 

(a;^ + gg)do(fe) + K(fe)-d^(fc)K(fc) 

[uj^ + Elik)] [uj'^ + E^{k)] 

(a;^+ad(fc)-[dg(fc)-d^(fc)]d*(fc) 

[io^ + Elik)] [uj^ + El{k)] 


(65) 

( 66 ) 

(67) 

( 68 ) 








with 


E±{k) = ^£l + Do{k)±D{k). (69) 

The real functions D, and D arise from the expres¬ 
sion 


A(fc)A+(/c) = i:)o(fc)i + D(fc)cr, (70) 

and have the following explicit forms: 

Do{k)=dl{k)do{k)+d*{k)d{k), (71) 

D(fc) = do{k)d*{k) + dl{k)d{k) + i[d(fc) x d*{k)], 
D{k) = ^/D^. 


normal state]: 


IC{q) = = ^ / idk)^[GMk)Gfi^{k - q) 


+ F^p(.k)FUk -q)- Gi°j(fc)G® (fc - q) 


(75) 


here the (0) superscript denotes the normal-metallic func¬ 
tions. 

At 9 = 0, the kernel gives the superfluid density ns- 
/C(0) = (76) 

TO 

Taking into account the explicit spin structure of the 
Green functions (|63ll and (l64l) . we find 


Note that relation F^p{k) = F^^{k) [see Eq. (l64)) ] is 
realized here since we explicitly consider the supercon¬ 
ducting state originating from Eqs. ([Ml) and ([M]), and 
having only principal components of the anomalous av¬ 
erages. 


E. Meissner kernel and superfluid density 


Assuming the London gauge, qA( 9 ) = 0, we find the 
linear response of the current to the external vector po¬ 
tential in the form 

Mq) =-^K:^j{q)Aj{q), (72) 

with the Meissner kernel 


F^ijiq) = — ({dk)^^dhl [Ga/3(fc)G/3a(fc - q) 
m J ml 


F Fap{k)Fp^{k — q) -\- (73) 


where q — (q, ti) with the bosonic Matsubara frequency, 

= 21 'kT, and n is the electronic density. 

The tensor structure of the kernel depends on the or¬ 
bital symmetry of the superconducting state. We dis¬ 
cuss here for brevity isotropic (s-wave) superconducting 
states, when the Green functions do not depend on the 
direction of the wave vector. Integrating over d^k in Eq. 
dZl, we choose q as the polar axis 2 of the spherical coor¬ 
dinate system. Then, due to integration in the azimuthal 
plane all nondiagonal components of /Cy (g) vanish. The 
tensor is thus diagonal, with K-xx = K.yy due to symme¬ 
try, while K-zz is generally different. At the same time, 
due to the London gauge, A( 9 ) does not have a compo¬ 
nent along q, so that the current (|72|) is insensitive to the 
K-zz component, which allows us to choose it arbitrarily. 
For simplicity, we take 


/Cy(9) =/C(9)%. (74) 

The integral in Eq. (1751) is divergent, and we employ 
the usual trick regularizing the divergency by subtract¬ 
ing the normal-metal expression [since lC{q) =0 in the 


n 


= ^E 


d^ 


Go(Ci ^) + Gl^(^, w) 


+ F*(e,a;)Fo(e,ca)+F*(e,ca)F(^,ca) . (77) 


Note the positive contribution from the anomalous com¬ 
ponents of the Green function, which stems from relation 
Fapik) = Fp^{k) [see Eq. (IMl) ]. which is realized here 
since we explicitly consider the superconducting state 
with only principal components of the anomalous aver¬ 
ages. 


F. Principal s-wave singlet even-ta state 


In the s-wave singlet (hence, even-w) case, the triplet 
component d(fc) is zero, and we have 

u^+e + \do\^' ° a;2 + e + |doP’ 

where dp = (io(|k|jw) is an even function of oj. Even 
though dp depends in general on jkj, we make the usual 
approximation of weak dependence near the Fermi sur¬ 
face, so that during integration over ^ in Eq. dZ?]) we can 
fix |k| = kp. Then, we reproduce the standard result 

Mop 


- = -Tj2- 3/2, 

n ^ {uj^ + \dp\^f^ 

and the superfluid density is positive. 


(79) 


G. Principal s-wave triplet odd-ta—dia state 


In the s-wave triplet (hence, odd-ta) case, do = 0 while 
d = d(|k|,a;) is an odd function of ui, and we have 


_ (iuj + -l- -l- Dp) 

° “ (a;2 +^2^Dp + D) {uj^ +e + Do-D)^ 

{iuj + 


Go = - 
G = 


(ca2 .^^■2^Dp + D) (u;2 + ^2 + Dg - D) ’ 

= 0 , 

F =_ 0^^ + + Dq _ 

(ta2 +^-i^Dp + D) (a;2 + Dp - D) ’ 


(80) 

(81) 

(82) 

(83) 
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where 


Do = dM, D = i [d X d*], D= yj - (dd)*(dd). 

(84) 

We again approximate d(k) by its value at the Fermi 
surface. Then, Eq. (EZl) yields 


n 


UJ 


1 


2Douj^ + 2Dl - 2D^ + DqD 
(w2 + Do - 


2DoUj^ + 2Dl - 2D^ - DqD 
(a;2 + Do + 


(85) 


Taking into account that D ^ Do [see Eq. il] , one can 
check that the expression in the square brackets is always 
non-negative, so ns > 0. In the case of unitary pairing, 
i.e., D = 0, the result has a form similar to the singlet 
case (17^ : 


ns 

n 


(jJ 


d*d 

(a;2 + d*d)^/^ ■ 


( 86 ) 


Thus, the superfluid density is positive also for the s-wave 
triplet odd-w case. 


H. Induced s-wave triplet odd-ca—para state 

The odd-w state discussed above [Eqs. (1501) - (1551) ] 
demonstrates diamagnetic response to the external mag¬ 
netic field {ns > 0), and, according to Refs. [13,11^, is 
therefore consistent as a principal superconducting state 
(while the question of finding a proper interaction is to 
date unclear). Below, we discuss the induced s-wave 
triplet odd-ca-para state (which is undoubtedly real ized , 
e.g., in S/F H m, m El or TS/DN [M El, [H [IMl 
proximity structures) using the same language. 

Since the odd-ca-para state is induced, it does not have 
a pair potential in the corresponding symmetry channel, 
while the corresponding superconducting correlations are 
described by anomalous averages. In order to use the 
language established above, we can still use notations A 
and A+, however, consider them now simply as auxiliary 
quantities parametrizing the Green functions. The Green 
functions found in microscopic models are, of course, 
mo del-dependent, but here we are interested in funda¬ 
mental properties determined by symmetries. For an ex¬ 
ample of a simple microscopic model, see Appendix [Bj 

As has been previously discussed [o^, , the odd-w- 

dia and -para states are characterized by different signs 
in the relation between F and D+. Instead of Eq. (1571) 
with SA = — 1 (odd-w-dia state), one has sa = 1 in this 
relation for the odd-w-para state (this state is realized 
in microscopic models with even-w pair potentials). In 
order to capture this property, we parametrize the Green 
functions taking the same form (EOl) for Aas{k), and de¬ 
fine A+ according to 

A+^(I,2) = -A^JI,2), A+^{k) = (87) 


note that the signs are different from Eqs. (EU and (l59l) 
(which were valid, in particular, for a principal odd-w- 
dia state). The Green functions can then be found from 
Eqs. (1531) and dsni)- The result can be obtained from Eqs. 
(I80ll - (l83l) by inverting the signs in front of Dq and D in 
all expressions. 

Instead of Eq. (I55|l we now have 


ns 

n 


^tY. 

to 


I 

Jd 


-2DoOj^ + 2Dl - 2D2 - DoD 
(a;2 - Do - 


-2DoUJ^ + 2Dl - 2D^ -b DqD 
(w2 - Do + Df/^ 


( 88 ) 


We assume that the w dependence of Do and D is such 
that w > VDo — D, so that the combinations in the de¬ 
nominators of Eq. (l88l) are positive (it can be verified 
that this phenomenological assumption corresponds to 
existing microscopic models for the odd-w-para state). 

Taking into account that D ^ Do [see Eq. (15^ ]. one 
obtains that the expression in the square brackets in Eq. 
(1551) is always non-positive, so ns < 0 (paramagnetic 
response). In the case of unitary pairing, i.e., D = 0, we 
find 


ns 

n 


UJ 


d*d 

(w2 — d*d)^^^ 


(89) 


III. COEXISTENCE OF ODD-to-DIA AND 

odd-lj-para states 

A. Superfluid density 

Now, we consider a possibility of coexistence of odd-w- 
dia and -para states in some region of space. We assume 
that the Green functions are linear combinations of the 
two contributions: 


GdijS — 

(90) 


(91) 

= (^/)«/5 + iF^)o.F 

(92) 


where the dia and para contributions, in accordance with 
our previous consideration, have the following properties: 

{F+U{k) = (D,)^„(fc), (93) 

{F^Uik) =-iF,);^{k). (94) 

Instead of Eq. dZZl), from Eqs. (173 and dZSl) we now find 


— — T I {God + Gop)^ + (Gd -|- Gp)^ 


+ {Fod + Fop){F*^ - F*^) + (F, + Fp)(rS - F;) . 

(95) 
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This result contains separate contributions from the dia 
and para states, as well as the cross term 


Then, instead of Eq. (1521) we can write the fermion part 
of the action as 


— = T ^ [2GodGop + 2GdGp 

n ^ J-oo 

+ (F*,Fop - FodFo;) + (F^Fp - F^F^)]. (96) 

We see that the contribution from the anomalous func¬ 
tions is purely imaginary, so the cross term is complex¬ 
valued. 

Considering unitary pairing for simplicity, we have 


God = - 


Frf = 


ioj 


-1- ^2 -1- 


Gop — -- 


iw -I- ^ 


(^2 -I- ^2 _ ’ 

(97) 

Gd = Gp = 0, Fod = Fop = 0, (98) 

dd ^ 


^2 -1- ^2 _j_ 


F — 

± rt — 


W2 -P ^2 _ 


Op 


-, (99) 


and then 


^ = i— 

n y Dod + Dop 


-|- 2Fod + ddd* — d^dp 


\/w2 -I- Dod 


4u!^ — 2 Dqp -|- ddd* — d^dp 


v'w2 - D, 


Op 


. ( 100 ) 


This expression is complex due to the purely imaginary 
combination (ddd* — d^dp). 

A complex-valued cross term would mean a complex 
current and is therefore unphysical. We conclude that 
assuming a possible coexistence of odd-w-dia and odd- 
Ld-para states, we arrive at an unphysical result. 


B. Josephson junction 

Now, we consider a Josephson junction in the tunnel¬ 
ing limit, assuming that the banks are represented by 
odd-o; superconducting states. Our main interest is the 
combination of dia and para states. 

The tunneling contribution to the action has the stan¬ 
dard form. 


So + 5aux + St 



ThRaPy 

Mnajd) 


II 

( "lrir2 
1 0 

II 

1 0 


0 

_ 

^rir2 

-r ) 


i5a/3(5(ri 

PH 

dapdiji 

PH 


7-2), 


r2). 


'f'/3(2), 

(103) 

(104) 

(105) 


After integrating over fermions we find 


So + ^aux + St 


— Tr In 

2 


(Ml 

\Xrl 


7lr\ 
Mr) ’ 


(106) 


where the matrix under the logarithm is written explic¬ 
itly in the LR space, while its elements are 4x4 matrices 
in the PH ® spin space. 

In the tunneling limit we expand the logarithm over 
the nondiagonal (in the LR space) part. The tunneling 
contribution arises in the second order: 

St = ^ Tr (MZ^UrMr^Trl) ■ (107) 

We make the standard assumption that the tunneling 
matrix element in the momentum representation 7kk' 
does not depend on the momenta (we denote this value 
by r= |r|e“). Then, 


= {d^kT){d^kR) 


xTr 


irP {Gl{uJ, kL)GRiuj, kR) + G'r{cO, kL)G'R{cO, kfl)} 


- r*"Fi(cd, ki)F+(cd, kR) - T^F+{uj, ki)F«(cd, k^) 


(108) 


where the Green functions are written as matrices in the 
spin space. 

For a junction between two conventional s-wave singlet 
even-w superconductors (with superconducting phases 
ifiL and (Pr)j we obtain: 


St 


lV'2a(l’L,T)V’fla(rfl,T) 


+ V‘^^^'llj*R^{rR,T)'ljjLairL,T)], ( 101 ) 


with the tunneling matrix element Ttltr- 

In addition to the particle-hole (PH) and spin spaces, 
we now also have the left-right (LR) space, so we deal 
with the direct product of the three spaces: PH 0 spin 0 
LR. In order to write the quadratic (over fermions) part 
of the action in a compact form, we define the vector field 


= (102) 


St 


ttG + IdpLlMofll cos 1/3 

2 ^ A/a;2 -I- |doLp\/w2 -|- |doi?p’ 


(109) 


where tp = {pR — pl + 2a) and G is the interface con¬ 
ductance in units of e^ jh, 

G = M\r\^VLVR, (110) 


with vl(r) being the normal-metallic density of states 
in the L(R) superconductor. The phase a of the tunnel¬ 
ing matrix element only shifts the superconducting phase 
difference, so we can set a = 0 and then deal with the 
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junction characterized by zero phase difference in equi¬ 
librium. The anomalous part of the tunneling action de¬ 
scribes Josephson coupling and can be written in terms 
of the critical current Ic'. 

5j = -^cos(p, (111) 

leading to the standard Josephson relation I = JcSint/?. 

Considering a junction between two triplet unitary (for 
simplicity) odd-w superconductors [see Eqs. 
we find 


{d^kL){d^kR) 


2\T\^GoL{ioM)GoRiuJ,kR) 

- SRT*^FL{uJ,kL)FR{^,kR) 

- SLT^Fliu;,kL)FR{uj,kR) 


( 112 ) 


where sr and sr are the signs originating from Eqs. dMl) 
and (IM)) and defined as follows: sr(^r'^ = ±1 if the cor¬ 
responding superconductor is odd-w-dia/para. Conse¬ 
quently, if both sides of the junction are of the same 
type, the Josephson (anomalous) part of the action is 
real, while if they are of different types, the Josephson 
contribution is purely imaginary. 

The d vector for the unitary pairing can be represented 
as with a real vector d. Then, taking the Green 
functions given by Eqs. dlzli-dMl), we obtain 

dRdRjsLe^^^ + SRe-^‘^)/2 

a;2 -h sz,d2 y/^2 




The phase-dependent combination in the anomalous part 
takes a form depending on the types of superconductors 
composing the junction: 


-I- src 
2 


— cos (/?, odd-w-para/para, 
cosifi, odd-w-dia/dia, 

—isin(/3, odd-w-para/dia. 


(114) 

We can obtain equivalent results in the quasiclassical 
technique (£-integrat ed Green functions) with proper 
boundary conditions [1021 . llO.jj . 

Strikingly, the last case (odd-w-para/odd-ta-dia junc¬ 
tion) yields purely imaginary Josephson coupling, leading 
to the Josephson current proportional to i cos ip. This 
result is unphysical and signifies problems regarding pos¬ 
sible coexistence of odd-frequency pairings with different 
types of Meissner response. 

Summarizing the results of the present subsection, we 
see that while a conventional even-w/even-w junction 
demonstrates the standard Josephson current propor¬ 
tional to sin:/:, the situation with odd-w superconduc¬ 
tors is more complicated. In even-w/odd-w junctions, 
the first-order Josephson coupling is absent due to differ¬ 
ent spin structure of the banks. Same-type combinations 


odd-w-dia/odd-w-dia and odd-ta-para/odd-w-para lead 
to sinusoidal current-phase relation (the general sign can 
be negative, which corresponds to the tt junction). At the 
same time, considering a different-type odd-ta-para/odd- 
ta-dia junction, we arrive at the Josephson current pro¬ 
portional to i cos if, which is imaginary and thus unphysi¬ 
cal. Similarly to Sec. lIII Al we conclude that the assump¬ 
tion of coexistence of an odd-ta-dia and an odd-ta-para 
state leads to unphysical results. 

One could try to avoid the contradiction, assuming 
that in odd-ta-para/dia junctions only configurations 
with di T d/j are realized, so that the imaginary contri¬ 
bution vanishes in the first-order Josephson coupling in 
Eq. (11131) . However, an argument of this sort (pair poten¬ 
tials adjusting to each other in a junction) could work if 
other configurations were energetically unfavorable (i.e., 
they would yield a free energy not in the minimum). In 
our case, they yield a complex free energy, hence are 
simply unphysical. Therefore, this argument can hardly 
solve the issue. 

Calculations for a specific microscopic model of the 
odd-w-para state in Appendix |E] confirm the results of 
the present subsection. 


IV. DISCUSSION 

Recently, it was demonstrated that in unifo rm sys¬ 
tems with broken time-reversal symmetry |104| and in 
nonuni form systems with preserved time-reversal sym¬ 
metry [lOSj j. even-w and odd-w superconducting compo¬ 
nents can mix with each other. At the same time, it 
turns o ut that th e states are divided into two separate 
classes fiol flO^ : while the even-w-dia state can mix 
only with the odd-w-para state, the even-w-para state 
can mix only with the odd-w-dia state. The two classes 
seem to avoid coexistence between themselves. 

Erom this point of view, when studying odd-w- 
para/odd-w-dia mixtures, we assume coexistence be¬ 
tween representatives of the two different classes. Our 
results demonstrate that this leads to unphysical super¬ 
conducting transport properties of the systems. 

How can this contradiction be resolved? Supercon¬ 
ductivity belonging to the first class is realized, e.g., in 
conventional s-wave superconductors (even-w-dia) and 
in proximity systems with conventional superconductors 
(even-w-dia and odd-w-para states), and these states are 
described by well-established microscopic models. At the 
same time, the question of possible realization of super¬ 
conductivity belonging to the second class is still open. 
The contradiction arising from the assumption of coex¬ 
istence between the two classes then raises doubts about 
the actual possibility to realize the second-class super¬ 
conductivity (in particular, the odd-w-dia state). 

On the other hand, if the odd-w-dia state cannot be 
realized, then what is possibly wrong in the arguments 
of Refs. [13, [H, where existence of this state was pro¬ 
posed from the general viewpoint of symmetry and sta- 
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bility (without presenting explicit microscopic interac¬ 
tion leading to this state)? The authors of Refs, [o^, 
argue that the odd-w-dia state can be realized only in 
a system with strongly retarded interaction, and can¬ 
not be described by a mean-field Hamiltonian. Assum¬ 
ing existence of a mean-field Hamiltonian immediately 
leads to EccJSSwith = 1 that signifies the odd-w- 
para state [^, [9^ . On the other hand, an effective re¬ 
tarded low-energy theory in the path-integral formulation 
[Eqs. (O-dl]) or (l33ll - (l^ ] emerges after integrating out 
high-energy degrees of freedom in a many-body Hamilto¬ 
nian. Assuming some general initial many-body Hamil¬ 
tonian H (that, e.g., contains all electrons and nuclei of 
the solid with their mutual interactions), one can define 
the Heisenberg operators ^/;Q,(r,r) = e^'^^/;Q,(i‘)e“'^’’ and 
'0+(r,r) = e^'^'0)j(r)e“'^’', and the anomalous Green 
functions 

Fa/3{ri,Ti-r2,T2) = {Tripa{ri,Ti)tl!i3{r2,T2)), (115a) 

^a/3(l‘l,n;r2,r2) = {Tri>^{ri,Ti)ljj^{r2,T2)), (115b) 

where averaging is over the exact many-body state. We 
are interested in the general relation between F and F '^, 
which can be obtained directly from definitions (11151) . fol¬ 
lowing the same logic as discussed for the mean-field case 
in Refs, Applying this to the homogeneous case, 

we immediately find the relation F^^{k, oj) = F^^ik, —uj) 
between the many-body anomalous averages, which has 
exactly the same form as Eq. (l57ll with SA = 1- This 
seems to be an unavoidable fundamental relation. There¬ 
fore, in order to obtain the odd-w-dia state described by 
relation (1571) with sa = one has to assume the im¬ 
possibility of a Hamiltonian description of the system at 
any level (both mean-field and many-body). This seems 
unnatural to us and increases our skepticism about a pos¬ 
sible realization of the odd-w-dia state. 

Having established a fundamental incompatibility be¬ 
tween a Hamiltonian description of the system at any 
level and the possibility of the odd-w-dia state, we still 
have to explain how relation (113 with Sa = can be 
avoided if one starts from the path-integral formulation of 
Eqs. (155)) - (1551) . We suppose that the key issue here is the 
spontaneous symmetry breaking. Taking a trial path for 
the pair potential instead of integration over the A* and 
A fields [introduced for the Hubbard-Stratonovich trans¬ 
formation in Eq. mi we arrive at the mean-field defini¬ 
tions of the anomalous averages, Eqs. m- This step in 
the derivation assumes that while the free energy has a 
superconducting manifold of equivalent minima (with ar¬ 
bitrary superconducting phase), due to an infinitesimally 
small perturbation (which is not even explicitly consid¬ 
ered) the system chooses some definite phase (thus break¬ 
ing the gauge symmetry). This choice of a single point 
from the manifold corresponds to taking the mean-field 
value of A (the trial path) instead of integrating over this 
field. 

At the same time, the symmetry-breaking (phase¬ 
fixing) perturbation is an essential issue. In the 


case of conventional superconductivity, we can write 
it explicitly in the second-quantized representation as 
/ dr[(5o'0|(r)V)j(r) -F (5^'0|(r)^^(r)], with |5o| 0. This 

Hermitian term breaks the global gauge symmetry, set¬ 
ting the preferential value of the phase to be equal to 
the phase of dp. On the other hand, for the odd-w-dia 
solution of Refs. [13, [Eq. (1571) with sa = —1], our 
hypothesis is that there is no spontaneous physical per¬ 
turbation (fluctuation) that can lead to the spontaneous 
gauge-symmetry breaking, thus fixing the phase of this 
state (this point of view correlates with the absence of a 
Hamiltonian description of the odd-w-dia state). Then, 
one has to retain integration over the superconducting 
phase, and the anomalous averages (FT)) vanish, so that 
Eq. (l57l) with sa = — 1 is trivially satisfied. In other 
words, the corresponding odd-w-dia minimum manifold 
can exist, but the symmetry cannot be spontaneously 
broken since there exists no physical fluctuation that 
could fix the phase. Note that in a related problem of 
an odd-w-dia state in Ref. [l^, the authors admit that 
the properties of the anomalous averages “depend on the 
properties of the unphysical, time-dependent, symmetry¬ 
breaking field.” 


V. CONCLUSIONS 

We have considered physical consequences of the 
odd-w-dia superconducting state proposed recently in 
Refs. [13,11^. This is the odd-frequency state with dia¬ 
magnetic Meissner response, which has no stability issues 
in the bulk and could therefore be realized as a prin¬ 
cipal superconducting state (at the same time, a con¬ 
vincing demonstration of such a state in a microscopic 
model is still lacking). Assuming the possibility to re¬ 
alize this state (due to a proper retarded interaction), 
we have studied its coexistence with the odd-w-para su¬ 
perconducting state, which is known to be generated 
as an induced superconducting component in, e.g., S/E 
or TS/DN proximity structures. Calculating the super¬ 
fluid density of the mixed odd-w-para/dia state and the 
Josephson current in the odd-w-para/dia junction, we 
find that the currents in both cases have imaginary con¬ 
tributions and are therefore unphysical. 

Taking into account rigorous microscopic derivations 
of the odd-w-para state in a number of models, we thus 
encounter the question of an actual realizability of the 
odd-w-dia state. Further analysis shows that a realiza¬ 
tion of this state implies the absence of a Hamiltonian 
description for the system at any (mean-field or many- 
body) level. Technically, the essential difference between 
the two states is described by different signs in the rela¬ 
tion between the two anomalous averages, see Eq. |57D 
that has sa = 1 for the odd-w-para and sa = — 1 for 
the odd-w-dia state. We conclude that in the latter case 
there is no physical perturbation leading to the global 
phase-symmetry breaking, and in the absence of a phys¬ 
ical perturbation in a mean-field Lagrangian formalism. 
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Eq. dSll) with sa = ~1 is trivially satisfied because the 
anomalous averages vanish due to integration over the 
superconducting phase. This also seems to be the only 
way to avoid unphysical imaginary current components 
in such a hypothetical odd-w-dia superconductor when 
brought in contact with the well-established induced odd- 
w-para state in, e.g., S/F or TS/DN hybrid structures. 
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Appendix A: Fourier transformation 

The Fourier transformation is introduced according to 
the following rules: 

A(l,2) = J{dk) j (Al) 
A{k, k')= A(l, (A2) 

where 

fc=(k, w), fca: = kr — WT, (A3) 

(A4) 

I 


The fermionic Matsubara frequencies are w = 7rT(2n-|-l). 
According to these definitions, 

^ ^ ^(27r)3j(k - k') = (5 ■ 

(A5) 

So, for the Fourier transforms of the Green functions 
and the pair potentials we have 


G^p{k,k') = -G'0^{-k',-k) 


= -{i’aik)i’}{k'))^p, 

(A6) 

Fap{k,k') = {'ipa{k)il)p{-k'))^^, 

(A7) 

F+^{k,k') = {i,:{-k)i,;{k'))MF^ 

(AS) 

A+^{k,k') = A*p^ik',k). 

(A9) 


Appendix B: “Superconducting ferromagnet” as a 
microscopic model of the odd-a;—para state 

In Sec. IIIIBl we discuss Josephson junctions between 
odd-w-para and odd-w-dia states introducing the former 
in a phenomenological manner with the help of auxiliary 
quantities A and A+, in order to use the same language 
as for the odd-w-dia case. We require Eqs. (|87)) to be ful¬ 
filled, so that the resulting state reproduces the symme¬ 
try and properties of the odd-w-para state known from 
microscopic models. At the same time, since clear micro¬ 
scopic models of the odd-w-para state are available, in 
this Appendix we check our general conclusions taking a 
specific microscopic model. 

The simplest microscopic example of the odd-w-para 
state is realized in a “superconducting ferromagnet,” 
which is the conventional singlet superconductor (pair 
potential A) with homogeneous exchange field h. The 
exchange field leads to appearance of the (induced) s- 
wave triplet odd-w-para component, in addition to the 
conventional s-wave singlet even-w-dia one. 

Instead of Eq. (IM)) we now have 

S0= J bP: {dr + 0 V'a + - ^IV'l)] , (Bl) 

and instead of Eq. (|53]) we obtain 


Mc./3(1,2) = 


(5(1 - 2) [Sap{dr + C) + H^3)al3] lA((T2)a/3 

-iA*{a2)aP ^(I - 2) [5al3{dT - 0 - h{^3)al3]) ’ 


The Green functions (matrices in the spin space) are then given by 

[{ioj + ^){oj^ + + |Ap — K^) + 2iwh^\ — h -I- J- |Ap — h'^ + 2iw{iw + ^)] da 


G{k) = -^ 


(a;2-he2-b|A|2-/i2)2 + (2a;/i)2 

F{k) = Fs{k) + Ft{k), F+{k) = F^Ak) - Flik), 


r 


(B2) 


(B3) 

(B4) 


where the singlet and triplet parts of the anomalous 


Green function have the following form: 


Fs{k) 

Ft{k) 


(w2+C2 + |A|2-/l2)A 

( w 2 + ^2 + I A |2 _ / j 2)2 + ( 2^/,)2 H ^ 2 ), 

—2iujhA ^ 

(w2 + ^2 + I A|2 _ h^)2 + (2^/j) 2 (*^3fT2). 


(B5) 

(B6) 
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Note that Fs{k) and Ft{k) have even-w/odd-w symme¬ 
try, respectively, and Eqs. (IB41) are consistent with the 
general relation (1571) at sa = 1 (since the pair potential 
belongs to the even-w symmetry class). 

Now, we consider a Josephson junction taking the 
above superconducting ferromagnet as the left bank of 
our junction, with A = |A|e®‘^^. At the same time, we 
assume that the right bank of the junction is the prin¬ 
cipal s-wave triplet odd-w-dia state, determined by Eqs. 
(Unil-dlSl). For simplicity, we assume the unitary case 
with d = and a real vector d, which is an odd 

function of w. Then, 


-F =- 

1^2 -I- ^2 -)- jJ2 1^2 _|_ ^2 _|_ (J2 

F = F(ido-2), F+= F*(i<T0-2)^ 


(B7) 

(B8) 


Although the anomalous Green functions of the left 
bank, Eqs. (IB4I) . are mixtures of the even-w and odd-w 
components, the lowest-order Josephson coupling in the 
junction is provided only by the triplet component from 
the left bank, so that the Josephson coupling is effectively 


between the odd-w-para and -dia states. 

The Josephson action [the anomalous part of the tun¬ 
neling action (11081) ] takes the form 

„ .ttG 

Sj = cos 95 

xy;im( 1^1 ) . (B9) 

Both the factors under the sum in the right-hand side are 
odd functions of w, so that their product is even and the 
sum is nonzero. Note that the phase dependence in the 
case of this odd-cu-para/dia junctions is determined by 
cosy) [in contrast to the last line in Eq. (11141) ]. because 
the phase of the odd-w-para component in the super¬ 
conducting ferromagnet is shifted by 7r/2 with respect to 
the phase of A [note additional i in the numerator in Eq. 
(lB6l) 1. 

At the same time, the main qualitative result of Eq. 
(IB9I1 is that the general consideration of Sec. 1111 Bl is con¬ 
firmed and the Josephson action for the odd-w-para/dia 
junction turns out to be purely imaginary and thus un¬ 
physical. 
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